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Abstract. We introduce an alternative definition of the relative height h'^{x) of a 
one-dimensional fluctuating interface indexed by a continuously varying real paramater 
< K < 1. It interpolates between the height relative to the initial value {i.e. 
in a; = 0) when k = and the height relative to the spatially averaged height for 
K = 1. We compute exactly the distribution P'^{hjn, L) of the maximum hm of these 
relative heights for systems of finite size L and periodic boundary conditions. One 
finds that it takes the scaling form P'^{hm, L) — L~^^'^ f'^{hmL^^^'^) where the scaHng 
function f'^{x) interpolates between the Rayleigh distribution for k = and the Airy 
distribution for k = 1, the latter being the probability distribution of the area under a 
Brownian excursion over the unit interval. For arbitrary k, one finds that it is related 
to, albeit different from, the distribution of the area restricted to the interval [0, k] 
under a Brownian excursion over the unit interval. 
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1. Introduction 

While the study of extreme value statistics (EVS) is by now a longstanding issue in 
the fields of engineering [1], finance or environmental sciences [3], it was recognized 
only recently that EVS plays a crucial role in the theory of complex and disordered 
systems [4]. Since then, the understanding of EVS in spatially extended systems has 
been the subject of many recent theoretical studies in statistical physics. Indeed, 
although the EVS of a set of independent or weakly correlated random variables 
is well understood, thanks to the identification of a limited number of universality 
classes, much less is known for the case of strongly correlated random variables. In 
that respect, the study of fiuctuating elastic interfaces has recently attracted much 
attention [5l[6l[7l[8l[9l[l0l[IIl[l2l[l3l[l4]. Although simpler to study, these models, 
in low dimension, furnish an interesting example of strongly correlated variables where 
analytical progress can be made. 

Here we are interested in the extremal statistics of the relative height of a one- 
dimensional elastic line described, at equilibrium by the Gibbs-Boltzmann weight 
Pst oc exp (— 7i) with 



where H{x) is the height of the interface and L is the linear size of the system. Such 
models ([I]) have been extensively studied to describe various experimental situations 
such as fiuctuating step edges on crystals with attachment/detachment dynamics of 
adatoms [Mj. Here we also impose periodic boundary conditions (pbc) such that 
H{0) = H{L) and the process defined by H in Eq. ([T]) is thus a Brownian bridge, i.e. a 
Brownian motion constrained to start and end at the same point. The Hamiltonian Ti. is 
obviously invariant under a global shift H{x) H{x)+c and therefore a more physically 
relevant quantity is the relative height. Here relative means measured with respect to 
some (arbitrary) reference value. Up to now two different reference values have been 
considered : (i) the spatial average value as in [H (Tj [HI [9l [TOl [HI [12], [13] or (ii) the 
height at X = 0, H{0) (or at any other point in space because of pbc) which is a rather 
natural choice in the context of time series [H] (note that a recent work considered, 
for weakly correlated random variables, yet another choice where the reference is the 
minimum value of the height field on [0, L] [I5] which we will not discuss here). This 
yields two distinct definitions of the relative height 



Once a reference value is chosen, one defines hm = niaxo<x<L (or hm = niaxo<x<L h^) 
and we ask the question : what is the probability density f^'^{hm,L)l It is 
straightforward to see from Eq. ([T]) that 
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(zz) {h%xo)h%xo + x)) = xo(^l- , (4) 

so that the variables h^{x) or h'^i^x) are obviously strongly correlated : the computation 
of hm is thus non trivial. 

The distribution f^{hm,L), studied numerically in Ref. [6], was then computed 
analytically by Majumdar and Comtet [Zl [8] who showed, for the model described in 
Eq. ([T]) with pbc, that P^{hm,L) = L~^/'^ f x^^y{hmL~^/'^) and computed exactly fAiry{x). 
Its Laplace transform is given by [7] 

/ /Airy(x)e-^^dx = j9v^ V e-""P'''2-^'' , (5) 

n=l 

where a„'s are the magnitudes of the zeros of the standard Airy function Ai(2;) on the 
negative real axis |T7|. For example, ai = 2.3381 . . ., a2 = 4.0879 . . ., = 5.5205 . . . 
etc 1^7\. It is then possible to invert the Laplace transform [81 ITS] 

/Airy(^) = ^ E e-'-/^^'J'U{-5/6, 4/3, bjx') , (6) 

X 

n=l 

where 6„ = 2a^/27 and U{a,b,z) is the confluent hypergeometric function [I7j. This 
inversion is useful as it can then be evaluated numerically, and plotted as in Fig. [TJ In 
particular, for small x it behaves as [8] 

/Ai.,(x)~ A«9/V5e-2"f/27x.^ as x-0, (7) 

oi 



while at large x, fAuy{x) ~ {72^/6/^)x^e-^'' [I9]. 

Interestingly, the Airy distribution in Eq. ([6]) also describes the distribution of the 
area under a Brownian excursion over the unit interval. We recall that a Brownian 
excursion is a Brownian motion starting and ending at the same point and constrained 
to stay positive in-between, on the unit interval [HI [20]. Surprisingly, this distribution 
arises in various seemingly unrelated problems in computer science and graph theory 

On the other hand, T. W. Burkhardt et al. focused on f^{hm,L) [H] which is 
different from the Airy distribution. Indeed, it is easy to see that in that case the 
distribution of hm is just the distribution of the maximum of a Brownian bridge on the 
interval [0, L]. Therefore P^(hm,L) = L^^ f^iay^hmL^^) where 

/Ray (a;) = 4 a; exp {-2x^) , (8) 

which is known in the literature as the Rayleigh distribution. In the literature on random 
matrix theory, fRay{x) (which is plotted in Fig. [T|) is also known under the name of 
"Wigner surmise" as it describes the distribution of the spacing between consecutive 
energy levels in the Gaussian Orthogonal Ensemble (GOE). In Ref. [H] it was also 
pointed out that the distribution of hm = maxo<x<L /^(a;), for generic translationnaly 
invariant systems, is related to the distribution of near extreme events, which has 
recently attracted some attention in statistical physics [22] . 
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Figure 1. The dashed green Hne represents the Airy distribution function = /Airy: 
this is the pdf of the maximum measured relatively to the spatial average value of 
heights of the interface on the interval [0, 1]. The solid red line represents the Rayleigh 
distribution function = /Ray: this is the pdf of the maximum measured relatively 
to the first point of the interface on the interval [0, 1]. 



Although the definition of these two relative heights are a priori only slightly 
diff'erent (see Eq. ((21)), the distribution of their maximum in Eq. ((61) and in Eq. ([8]) has 
completely different expressions. These two definitions of the maximal relative heights 
thus lead to qualitatively different distributions : this is particularly striking if one 
considers their small argument behavior in Eq. (jTl) and Eq. (Ej). The goal of this paper 
is to understand how the distribution of the maximum evolves when one changes the 
definition of the relative height from h^{x) to h^{x). In that purpose, we choose an 
alternative reference value, indexed by a continuously varying parameter < k < 1 to 
define the relative height h^{x) as 

h^ix) = H{x) / H{x') dx' , (9) 

i^L Jq 

such that lim^^i = h^{x) and lim^.^oh'^^x) = h^{x) : h'^{x) thus interpolates 
between h^{x) and h^{x). Using path integral methods, we compute exactly the 
distribution P'^{hm,L) of the maximal relative height hm = ma.Xo<x<L h'^{x) and find 
that P'^{hm,L) = L~2f'^(^hmL~2) where f^i^x), which interpolates between the Airy 
distribution for k = 1 and the Rayleigh distribution for k = is given by 



-bn{q,K)/x^ 



xf/(-5/6,4/3,6„(g,/€)/x2), (10) 
where &„(g, k) = (2K/27)(a„ + (1 " i^)q^ / ^iid where the coefficients F„(g) are given 



MRH of one- dimensional interfaces : from Rayleigh to Airy distribution 



5 



below f l42i ). In particular, for small argument, we show that it behaves as 

The paper is organized as follows. In section II, we illustrate the method of path 
integral to compute the distribution /"(x), i.e. the distribution of the maximum of a 
Brownian bridge. In section III, we use this path-integral formalism to compute the 
distribution of the maximum of the relative heights h'^{x) defined in Eq. ((91). In section 
IV, we analyse in detail the limits k ^ 1 and k — of the scaling function /'"(x) in Eq. 
(fTOl) and in section V, we present the comparison of our exact results with numerical 
simulations for different values of < k < 1. In section VI, guided by the fact that 
/^(x) = /Airy(a^) describes the distribution of the area under a Brownian excursion over 
a unit interval, we show that /'^(x) is related to, albeit different from, the distribution 
of the area restricted to the interval [0, k] under a Brownian excursion over the unit 
interval. Finally, we conclude in Section VII. For clarity, some technical details have 
been left in Appendix. 



2. Distribution of the maximum of a brownian bridge : a path-integral 
approach 

In this section, we derive the distribution of the maximal relative height = 
maxo<2:<L /i°(x) where /i°(x) = H{x) — H{0). The joint distribution of these heights 
h^{x), given Eq. ([T]), is 

= exp [-\ £ ' dx j 6ih\0) - h^imh^O)) , 

(12) 

where is a normalization constant. In this expression (fT2ll . the first delta function 
ensures the pbc while the second one ensures simply the definition of h^{x) in Eq. ([2]). 
From this statistical weight, one directly obtains the cumulative distribution of the 
maximum F'^{hm, L) = Prob[maxo<a<L h^{x) < hm, L] as 

h°{L)=0 

F\h^,L) = ^^ I Vhix)e-yoi9^'^'(-)f X n 0{hm~h'{x)), (13) 

/iO(o)=0 0<x<L 

where product of Heaviside step functions constraint the paths to stay under hm. Note 
that the normalization constant Zi has just the same path-integral representation as in 
Eq. f |T3l ) but without the constraint ensured by the step-functions 

hO{L)=0 

Zl= J Vh{x)e-yoi9.hOi.)f _ (^4^ 

/iO(0)=0 



MRH of one- dimensional interfaces : from Rayleigh to Airy distribution 



6 



One sees explicitly on these expressions (fT3l fT4ll that F^{hm-, L) is the cumulative 
distribution of the maximum of a Brownian bridge on the interval [0, L]. One can 
now use path-integrals techniques to write F^{hm,L) as 

FV,„.L)^ (0|exp(-Lfl„)|0) 
(01 exp(-Li/„„)|0)' 

where -fffree = cind Hm = + ^m{x) where Vm{x) is a confining potential 

such that Vm{x) = if — oo < x < and Vm{x) = +oo if x > hm- The numerator 
and the denominator in Eq. (fTSll can be computed in a similar way: if one denotes 
\E) the eigenvectors of a Hamiltonian H, with eigenvalues E, H\E) = E\E) then the 
propagator can be expanded on the eigenbasis as 

{y\eM-LH)\x) = Y,rE{y)M^) e"^"" , (16) 

E 

where iPe{x) = {x\E). Applying this formula (fTBl) to Hm one has 

(0| exp(-L^^) |0) = -j^ dq{sm{q hm)fe-'^ = -j== (l - e~^^\l7) 

and similarly (0| exp(— L^free) |0) = (27rL)~2. And finally one obtains the cumulative 
distribution function as 



2 



F\hm,L) = l-e-'^ , (18) 

from which one deduces that the probability distribution function P^{hm,L) = 
dh^F^{hm, L) takes the scaling form (valid for all L for this model) 

P\hm, L) = , fix) = Axe-'^' , (19) 

which is known in the literature as the Rayleigh function. Of course there are other 
simpler methods to compute f^{hm,L) for the present case (see for instance [T4]). 
However the path integral method presented here can be extended to more complicated 
situations, and will be used below. 



3. Distribution of the Maximum Relative Height P'^{hm,L) 

In this section, we extend the calculation presented above to the computation of the 
cumulative distribution of the maximal relative height hm = maxQ<^^L h'^{x) where 
h'^{x) is defined in Eq. (I9|). As above we start with the joint distribution of the heights 
h'^{x), which from Eq. ([T]), can be computed as 



= — exp ( / ( ) dx I 5ih^i0) - h^iL)) 




'0 



h\x')dx' , (20) 



MRH of one- dimensional interfaces : from Rayleigh to Airy distribution 7 

where is a normalization constant which ensures that / Vh'^{x)P[{h'^}] = 1. The 
cumulative distribution of the maximum F'^{hm,L) can thus be written as 

h'^(L)=u 
J-oo J 



Jk •) -co 

/i«(0)=« 



Yl e{hm - h^'ix)) h^{x')dx''] , (21) 



X 

0<x<L 



which is the starting point of our exact computation. 
3.1. Normalization constant 

As above (see section [2]), the normalization constant has the same path integral 
expression as Eq. (l2Tl) . but without the constraint that each path must stay below hm, 
that is 

du / m"(a;)e-i/o^(^-'^''(-))'d-5 ( / h''{x')dx'\. (22) 

/iK(o)=u 

Let us make the change of variable y{x) = h'^{x) — u to obtain 

Z. = j du J Vy{x)e--^^o(d.y{-))''^^§ l^j^ y{x')dx' - kLu^ . (23) 

°° y{o)=o 

The path integral is, up to a normalization constant, the probability that the truncated 
area = J^^^ y{x)dx under a Brownian bridge {y{x), < x < L} is equal to kLu. 
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One thus has 



y{L)=o . 
J Vy{x)e~-2 fo(dM-))''i^S y{x')dx' - kLu 

Proba(ylK = kLm, L) = ^^^^ — 



y{L)=0 

J P?/(x)e"^-^o^(^^^(^))^'^^ 
y{o)=o 

The denominator is just (27rL)~^/^, as computed before (in Section[2]). Then we simply 
have 



hoc 



= —= I du Proba(ylK = kLu, L) (25) 



Changing v = kLu, and using the property that -P(Ak, L) is, by definition, normalized 



2411 , we find that 

^ 1 
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3.2. Cumulative distribution of the maximal relative height (MRH) 

Proceeding as in the first section, we now compute the cumulative distribution of the 
MRH using path integrals. We start with the expression given in Eq. (l2Ti l 

h"-(L)=u 
J — oo J 

X Yl e{hra - h%x))S ( [ dx'h%x')] , (27) 



where is given above (l26ll . The first integral in Eq. (l27ll indicates that we sum over 
all initial points u (which has to be below km), and then we sum over all paths beginning 
in u and constrained to stay entirely under the height km (imposed by the product of 
step functions). Performing the change of variable y{x) = km — h'^{x) and v = km — u, 
we get all /im— dependence in the delta function, and we constraint the paths y{x) to be 
positive: 

+00 y{L)=v 

F^{h^,L) = kLV2^ j dv jvy{x)e-^^^oMd.yf 

y{0)=v 

x6(k^KL- r\{x')dx'] H 9{y{x)). (28) 

\ Jo / o<x<L 

The constraint of positivity can be incorporated rather easily by adding a potential term 
Vo{y) in the exponential, with Vo{y) = +00 ii y < and Vo{y) = if ?/ > 0. Then we 
have 

+00 y(L)=v 

y{0)=v 

X 6 (^Lhjn — j y{x)dx^ . (29) 

Using 5{ax) = {l/a)6{x), we isolate the variable km, and we write 

+00 y(L)=v 

F^{hm,L) = VM j dv y"p2/(x)e-^o^Ml(^-?')'+^o(j/W)} 

y{0)=v 



X 5 - — y{x')dx'j , (30) 

and all the dependence on is contained in the 5-function. Thus its Laplace transform 
with respect to hmj F'^{p,L) = dhmF'^{hm,, L)e~P'^"' (notice that km is necessarily 
positive) has a simple expression 

+00 y(L)=v 

F^ip,L) = V2^ J dv Jvyix)e-^oH-.(9^^)'+My(^))}e~^fo"y(^^''^ . 

y(0)=v 
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(31) 

Let us call u = y{K,L) the intermediary point and using the Markov property, we can 
separate the path integral in two independent blocks: one for the time interval [0, kL] 
and the other for the time interval [kL, L]. We have 

+00 +00 y{KL)=u 

F^{p,L) = V2^ j Av j Au j py(x)e-^o''^'^45(9-!^)'+^o(y(x))+;^j/W} 

y(0)=v 

y{L)=v 
y(KL)=u 

The first factor is the propagator between time and kL of a quantum particle in the 
triangular potential Vo,p/(KL)(a;) = Vq{x) + ^x, denoted as G^iryiu, nL\v,Q) (the index 
'Airy' refers to the fact that the solution of the Schrodinger equation with a triangular 
potential Vo,p/(KL)(a^) is solved by Airy functions), and the second term, denoted as 
Gq{v^L\u^kL) is simply the propagator of a quantum particle in the potential Vq{x) 
between the times nL and L. One thus has 

F''{p,L) = V2^ j J {Go{v, L\u, kL) X GAiry{u,KL\v,0)} , (33) 



where the p— dependence is entirely contained in the Airy propagator. Note that our 
calculation is completely justified when < /t < 1, because of the cut in two parts of 
the path integral. We will examine in more detail the two limiting cases k — and 
K — > 1 below. 

The probability density function of the MRH P'^{hm,L) is simply the 
derivative of F'^{hm,L) with respect to hm- Its Laplace transform P^{j),L) = 
e-ph^p-{hm, L) dhm is thus P^{p, L) = pF-{p, L), 



hoo +00 



P^{p,L) =pV2^ j j (^u{Go{v,L\u,kL) X GAiryiu,KL\v,0)} .{34) 



We first start with Go{v, L\u, kL) which is the propagator of a particle with a wall at 
the origin. It is given by 

Go{v,L\u,K,L) = e - e (35) 



^27rL(l - k) 

2 /'+°° L{1~k) 2 

= — / dg sin(gt>) sin(g'u) e 2 , (36) 
Jo 

where the expression in the second line will be useful in the following. The 
propagator G'Airy(^, ^^-Z^l^^, 0) can be computed using the formula in Eq. (fT6l) (see also 
in Appendix A. 4 , formula (IA.25I1 ): 



MRH of one- dimensional interfaces : from Rayleigh to Airy distribution 



10 



GAivy{u, kL\v, 0) 

,-a„2-i/3(pv^)2/» 



UiJ h (Ai'[-a„l)' 



V — Or 



X e 



(37) 



We then perform the changes of variables v (2p/(/tL)) ^^^v and u (2p/(/tL)) ^^^u 
followed by g — {2p/ {KL)Y^^q to obtain 




oo oo oo 



P'^{p, L) = p\l — J du J dv J dg sin(gM) sin(gf ) exp -2 ^/^(l - /t)/? '^^^{pVLf/'^q^ 





n=l 



Ai(u — Q;„)Ai(w — an) 



[Ai'(-a„)]^ 



exp 



-a„2-i/3/s:i/3(p^)2/3 



(38) 



One recognizes the scaling variable pyL in the Laplace space which implies that the 
scaling variable in the real space is, as expected, hm/ V^, and then we can write the 
scaling law 



P\h^,L) = —f'^ 



hr 



where the Laplace transform of f^i^x) with respect to x is given by 



(39) 



oo oo oo 



/"(p) =p\ — I dq f du I dv sin(gu) sin(gt>) exp 



2^ ' V2 





oo 



X 



^ Ai(^-a„)Ai(.-a., [_«„2-V3,iV/3] _ 



(Ai'[-«„])^ 

This formula can be written in the more compact form 



(40) 



~np) 



Jdgf; [F„(g)]Ve-^"^''^'^'("-^''+(^-'^)'^"'^''''), 

n n=l 



(41) 



with 



Fn{q)= / dw^J^JJ — ^sin(gM). 
Ai -a„ 



(42) 



To invert the Laplace transform in Eq. (|4T1) we will make use of the result found by 
Takacs in Ref. [T8l 



27Tpe 



_2-l/3p2/3^^ 



(X) = ^e-''"/^%l/'U{-5/6,4/3,bJx^), (43) 



a; 



10/3 



where 6„ = 27^/27 and U{a,b,z) is the confluent hypergeometric function [17j. This 
yields here 



■ X 



"10/3 



TT 



dqJ2[U<l)? hn{q,Kf''e- 



bn{q,K)/x'^ 



n=l 



X f/(-5/6,4/3,6„(g,ft:)/x2). 



(44) 
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where bn{q,i^) = |f (o^n + ^^Q'^)^- This class of distributions indexed by the real 
parameter < k < 1 smoothly interpolates between the Rayleigh distribution for 
K — > and the Airy distribution when k — 1. This will be demonstrated analytically 
(in particular the behavior when k ^ is not obvious on Eq. (I44l) ) in the next section 
and numerically in section VI. 



4. Recovering Rayleigh and Airy distributions 

In this section, we show how f^{x) yield back - as it should - the Rayleigh distribution 
(respectively the Airy distribution) in the limit k ^ (respectively in the limit k ^ 1). 



4.1. The limit k ^ and the Rayleigh distribution 

We first study the limit k 0. From the definition of h^{x) given in Eq. ((91) one 
expects that it yields back the definition of h^{x), i.e. the relative height relative to 
H{0) in Eq. (I2|). This is also rather obvious on the joint distribution in Eq. f l20l) . together 
with oc 1/kL, that it yields the joint distribution of a Brownian Bridge. Therefore 
one expects indeed that /'^(x) will converge to the Rayleigh distribution /Ray(a^) in Eq. 
(E]). However, this limit k — in our formulae for f'^{x), or in its Laplace transform 
f^ip), is not trivial. This is partially due to the fact that the computation of Rayleigh 
distribution has been performed in "direct" space while for the Airy distribution, the 
computation is more conveniently performed in Laplace variable. Indeed, one needs an 
identity about the behavior of the Airy propagator when the time difference tends to 
zero. In [Appendix A.5[ we show the identity (1A.37I 1 which we then use straightforwardly 
in Eq. (I40ll to compute the limit k ^ 0. 

Indeed, from the Laplace transform of the scaling function in Eq. ( l40ll . one performs 



the change of variables v' 
we have 



(2p/k) 



-1/3 



f , u 



{2p/k) ^/^m, and q' = (2p/fi:)^/^g, so that 



rip) 



PM- 



X 



dq' I du I di;' sin(g'u') sin 



1/3 00 Ai 

E- 

n=l 



U 



a. 



Ai 



exp 

2p\ 1/3 



■2(1 - 



(45) 



(?) 



Or 



,2-1/3^1/3^2/3 



(Ai'[-«„])^ 

Now, the second line in that equation has a good limit when k — > 0, as shown in 
appendix [Appendix A.5[ We then make use of the identity in Eq. ( 1A.37I ) to obtain the 
limit K — ^ (with u = u' , v = v' , q = q') 



lim/"(p) 
«— >o 



dq 



du 



P\l- 



du e 



-pu 



dv sin(gM) sin(gt>) exp 
dq sin^(gM)e^'^ = 



q' 



X 6{u — v)e~ 



pu 



du Aue 



-pu 



(46) 
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Now we recognize immediately the scaling function in real space, since = 
/p°°du/°(u)e-P" such that: 



lim f'ix) = fix) = 4x6-''' , (47) 

K— >0 



as expected. 



4-2. The limit k. ^ 1 

The analysis of the limit k — ^> 1 is straightforward and can be done directly on the 
expression of our scaling function /"^(x) in Eq. (I44l l. In the limit k, ^ 1 one has 

lim bn{q; n) = bn{q; 1) = &«, (48) 

independently of q. Therefore the only dependence left in the integrand of Eq. (l44l) is 
contained in Fn{q). It is straightforward to obtain the useful relation 

oo oo oo oo 

|dg[F„(g)]^ = jdqjdujdv ^^^^ ^<<i^)^^^^ ^<<1-) 





oo 





where, we have used in the last line, the identity du{k\{u — q;„))^ = (Ai'(— a„))^. 
Using this relation (l49l) in the expression of f^ix) we get 

oo 

limr(x) = /Airy(x) = 2v^a;-i°/3^62/3g-f>„/..^ f/(-5/6, 4/3, fcjx^), (50) 

^ n=\ 

which is the Airy distribution function. 

5. Some properties of the distribution /"(x). 

In this section, we derive some properties of the distribution f^ix). 

5.1. Asymptotic limit of the scaling function 

In this subsection, we study the asymptotic limit of the scaling function when x <^ 1, 
which corresponds to h^ <^ L. One main difference between Rayleigh and Airy 
distributions is the singularity of the latter for x = 0. The Rayleigh distribution is 
perfectly analytic when x — (~ x), in contrast to the Airy distribution which exhibits 
an essential singularity /Airy ~ x~^e~^"i/^^^^, see Eq. (I7|). Here we show that /'^(x) 
exhibits an intermediate behavior with an essential singularity when x — > (see Eq. 
dnSl) below). 

We start with the expression for f'^{x) given in Eq. (I44ll which we recall here 

r(x) = ^x~^°/^ rdqf][Fn{q)f bn{q,^)'/'e-'-^'"^y^' 
f/(-5/6,4/3,6„(g,«:)/x2). (51) 
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that we want to analyse in the asymptotic limit x — 0. To this purpose, we perform a 
change of variable q = xq such that bn{xq, k)/x^ has the following small x expansion 



bn{xq, k) 



2Kar 



+ 



6ai{l 



K)q 



+ 0{x^ 



a;2 27x2 ' 27k 

so that, due to the exponential factor exp (— 6„(xg; k)/^^) in Eq. 
with n = 1 can be retained in sum over n. Besides, using that f/(— 5/6, 4/3, ~ y 
for 1/^1, [It] we have 



(52) 

5Ti ). only the term 

5/6 



r(x) 




with 



97r 



K 



1 



K 



3/2 



X 



exp 



K 



(53) 



ai)/Ai'(-ai; 



1.8173. 



27 x^ 

where we have used the small 
argument behavior of Fi{q) given in Eq. (l42l) Fi{xq) = x{c(i){^(l) + O^x"^). Of course, 
such formula is only valid for < k < 1 strictly. One sees on this expression f l53ll that 
this probability is exponentially small as soon as hm -C VkL. This behavior might be 
qualitatively understood by noticing that the configurations for which hm <^ VkL are 
actually essentially fiat on the interval [0, kL]. 

A comparison between this asymptotic behavior in Eq. (l53ll and the small argument 
behavior of the Rayleigh distribution (JHl) suggests that there exists a scaling regime 
X — s> 0, K — i> keeping the ratio x/y^ fixed such that 



f^ix) (XX g' 

with the scaling function 

constant 

9[u 



(54) 



u 



exp 



2al 1 



when u 
when u 



oo, 



0. 



(55) 



27 n2 

The limit m ^ oo means that, for x small but fixed, we take formally k — 0, and then 
we recover the asymptotic behavior of the Rayleigh distribution. On the other hand 
for M — > we keep x much smaller than k (actually x <^ -^/k as done previously), and 
one recovers the result above (1531) . Similarly, a comparison between the behavior in 
Eq. ^ and Eq. f l53ll suggests that there also exists a scaling regime where x — >^ and 
(1 — k) — > keeping x/ \^T^^ fixed such that 



/"(x) oc X ^ exp 



2af K 



h 



X 



where the scaling function h{v) has asymptotic behaviors 

constant when v ^ oo, 
when v ^ 0. 



h{v) 



(56) 



(57) 



Again w — > oo means that we keep x fixed, but small, and we take k — > 1, thus we 
recover the asymptotic behavior of the Airy distribution. When f ^ 0, we look at the 
region where x <C \/l — n <^ 1, thus recovering the behavior obtained above (l53l) . 
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5.2. Moments of the distribution. 

If one denotes the moments of /"(x) as = dxx'^f'^{x), it is straightforward to 
show that, for the Rayleigh distribution, M^^^ = 2~"/^r(l + n/2). On the other hand, 
the computation of the moments of the Airy distribution M^^^ is highly non trivial. In 
Ref. [IB], Takacs found a recursive method to compute them and recently, M.J. Kearney 
et al. found an explicit integral representation of these moments as [23] 

A /~ \ O poo 3{n-l). 

where Ai(2;) and Bi(z) are the two linearly independent solutions of Airy's differential 
equation y"{z) — zy{z) = 0. 

The computation of the moments for arbitrary n seems to be highly complicated 
and we have been able to compute only the first one = xf^i^x) dx. Indeed, 
permuting the integral over x and the integral over q together with the discrete sum 
over n and then performing the change of variable x x^Jhn{q, k,) one obtains that 
is actually independent of k, yielding 



Mr = ^y|, V«:G[0,1]. 



(59) 



6. Numerical results 
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Figure 2. Left : Test of the scaling relation with numerical plots. From the raw 
histogram P^^^{h,n, L), we draw as co-ordinate x = {hm + a)/VL, with a = 1/^/2 and 
as ordinate /'^(x) = VLPnuinihm., L). Histograms with different L but same k collapse 
in the same curve, which proves the scaling relation. Right : Test of the analytic result. 
The dots are the numerical simulations (where we use the 'translation' parameter a) 
and the curves represents numerical computations of the analytic function with the 
help of Mathematica. 



We have compared our analytical results with numerical simulations. To generate 
the configurations of the interface distributed according to the Gibbs-Boltzmann weight 
Pst oc exp (— 7i) with H given in Eq. ([T]), one could of course use standard Monte Carlo 
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method. However, this suffers from critical slowing down (with dynamical exponent 
z = 2) and this is not very efficient. Instead, as in Ref. [2lj, one can use the 
property that the process defined by this statistical weight Eq. ([1]) together with 
periodic boundary conditions is a Brownian bridge. To generate it, one first generates 
numerically an ordinary Brownian motion Bq = and Bi = -Bj_i + r]i where rji's 
are independent and identically distributed (i.i.d.) random variables drawn from a 
distribution Proba(?7 = x) = g{x). One can then generate a Brownian Bridge through 
the relation Hi = Bi — {ij L)Bl, which ensures pbc. One can show that this procedure 
yields the correct statistical weight for large system size L (when g{x) is a Gaussian, 
this procedure is exact for all L). The relative height is simply = Hi — Si=i 
and the distribution of hm = maxo<i<L is then computed for different system sizes 
L = 128, 256, 512 by averaging over 10^ samples. From the raw histograms representing 
Pnumi^m, L) , foT & givcu L, wc waut to cxtract the scaling function f^^^{x) to first check 
that the data for different sizes L, with k fixed, satisfy the scaling form that we have 
obtained analytically in Eq. ( l39l) . If we plot directly \fLP^^^{hm-i L) and as a function 
of a; = hm/ V^, one can observe some finite size corrections to this scaling form. 

These finite size effects can be understood using the analysis done in Ref. [I2] 
where the first finite size corrections to the Airy distribution were computed (see also 
the discussion in Ref. [25]). There it was found that 

^-"-("- ^' - iz (^) - ^/ (^) + " (^"'O 

where a is the non-universal constant (that depends on the details of the Hamiltonian), 
and = /Airy the Airy distribution function. This first correction can be resummed as 

PL.....{hm. L) = ^ (/^ + O (L-^)) , (61) 

which is a good trick to get rid of the leading finite size effects. Following this exact 
result for K = 1 ( !6Ti ). we thus plot on Fig. m (left) ^P^^^{hm, L) function of 
X = {hm + a)/ with a = 1/a/2 for all k and observe a very good collapse of curves for 
different L. In Fig. [2] (right), we show a comparison between \/LPI^^^{hm, L) and our 
analytical predictions in Eq. (I44ll evaluated with Mathematica for different values of 
K = 0, 1/32, 1/2 and n = 1. The very good agreement between analytics and numerics 
is obtained here via a single "fitting" parameter a = l/\/2 for all values of k, adjusted 
to take into account finite size corrections fl6TI) . 

7. Truncated area under brownian excursion 

In this last section we explore further the connections between the distribution of the 
MRH for arbitrary k and the distribution of areas below constrained Brownian motions 
[20] . Indeed, for k = 1, it was shown in Ref. [3,18], that f^{x) describes the distribution 
of the area under a Brownian excursion on a unit interval. Here we consider a Brownian 
excursion, i.e. a Brownian path {x{t), < r < T} starting and ending at the origin 
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kT 



Figure 3. Truncated area (in blue) under a realization of a Brownian excursion (in 
red). 



x(0) = x(r) = and which stays positive in-between Vr g]0,T[, x(r) > and consider 
the random variable = J^^ x{T)dT, which we call the truncated area (see Fig. [3|). 
We compute its distribution, which can be obtained by a similar calculation as shown 
above. 

We start with the joint distribution of positions of the brownian excursion 

P[x{t)] = MmJ—e-yoidrxf'-SixiO) - x{T))6{x{0) - e) ff 0(x(r)) 

= lim-^e-^o'M5(9--)'+vbWr))}^^^^Q^ _^^y^^^^^^Q^_^^ ^^g2) 
e^o Z[e) 

where e is a regularization parameter which is needed because of the non-smoothness of 
the Brownian path (see Ref. [ZlIH]) and Z{e) is the normalization constant 

x(T)=e 

Z{e)= j Px(r)e-^o''^Hl(«--)'+^oW-))} =Go(e,T|e,0) , (63) 

a:(0)=e 

where (^0(^5 to) is the propagator of the free particle with a wall in x = (l35l) . The 
distribution of the truncated area JJ^^ x(r)dr, for the Brownian excursion {x{t), < 
r < T} is then 

P( A, T) = Prob (^j^^ x(r)dr = A^j = (d (^A^ - j^^ x(r)dr^ ^ , (64) 
where the mean value is taken respectively to the distribution of excursions P[x{t)]. We 
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have 



J Px(r)5(A.-/Jx(r)dr) e-/o"dHl(9-)^+v;,(.M)} 
P(A,T) = lim^^^^^^ -— .(65) 

x{0)=e 

We now take a Laplace transform with respect to Ak, which is a positive quantity. So 
we have 



POO 

Pip,T)= / dAP(A,T)e-^-^'' 
Jo 



(66) 



x{T)=e 



J Vx{t) e-/o'dr{i(9.x)2 + Vb(x(r))} j-f Arx{r) 



By similar techniques, one can identify a convolution of propagators. Expanding 
numerator and denominator in powers of the regulator e, one finds that the lowest 
order is for the numerator and the denominator. Taking the limit e — 0, one finds 
that the distribution of the truncated area obeys the following scaling form 

m..r) = pip-«(^). (68) 



with the scaling function, parameterized by k, given by 

4^6 _io/3 fbniq,K 



^.^(a;, = a; -10/3/ dg^[gF„(g)] /^^^i^^^ .-b.M/i^^^) 



n=l 



2/3 



f/(-5/6,4/3,6„(g,ft:)/(/€x2)), (69) 

with the same coefficient bn{q,i^) as in Eq. (I44ll and -F„(g) defined in Eq[42l. We can 
check easily that the scaling function w'^{x) converges to the Airy distribution function 
as K ^ 1 , as shown in Refs [3, [8]. However, for k 7^ 1 these two distributions f^^^^ and 
w'^{x) are different, albeit having similar expressions. 



8. Conclusion 



To conclude, we have introduced an alternative definition of the relative height h'^{x) 
of an elastic one-dimensional interface with pbc, indexed by a real n G [0,1], which 
interpolates smoothly between the height relative to the initial value when n ^ and 
the height relative to the spatial average value for k ^ 1. We have obtained, using path- 
integral techniques, the exact distribution of the maximal relative height f^i^x) which 
interpolates between the well known Rayleigh (k = 0) and Airy (k = 1) distribution (see 
Fig. [1]). This thus constitutes one new family, parameterized by k, of strongly correlated 
variables where such analytical calculations can be done. Although our calculations have 
been done for a continuum model ([T]), the arguments presented in Ref. [T2], based on the 
Cental Limit Theorem, show that our results actually hold in the limit of large system 
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size, for a wide class of lattice models of interfaces H = J2i \Hi — -f^i+i |^ with arbitrary p 
and in that sense this distribution /''(x) in Eq. (I44l) is universal. Finally we have shown 
that the method employed here to compute this distribution can be used to compute 
the distribution of the truncated area (see Fig. [3l) under a Brownian excursion. 



Appendix A. Quantum mechanics problems 

In our path integral calculations of distribution of MRH, we have to deal with simple 
one dimensional quantum mechanics problem. This appendix summarizes the results 
used in the main text. 



Appendix A. 1. Free particle 



Let us consider a free particle in one dimension. The Hamiltonian is: 

1 d^ 

The solution of the Schrodinger equation are the planes waves: 

M^) = e^'^ (A.2) 

with — oo < k < +00. The associated energy to ipk is = The calculation of the 

propagator is done by decomposing onto the eigenbasis: 

G'f,ee(2/,t|x,s) = (y|e-(*-^)^f-°|a;) 



+0O 



VMt-s) V 2(t-s) 
Another way of doing it is to follow Feynman's prescription [26]: if the lagrangian is a 
quadratic form of its variables, then the propagator is proportional to the exponential 
of minus the classical action joining the two points in space-time. Moreover, if the 
lagrangian does not depend explicitly on time, the constant of proportionality is simply 
a function of the difference of the two times: 

G{y, t\x, s) = F{t - s) exp (-^[sciassicaKT)]) . (A. 4) 

When we are in imaginary time, because we are dealing directly with probabilities, the 
function F{t — s) can be obtained by normalization. For the example of the free particle, 
the action is: 

1 fdx(T)^ ^ 



S[xiT)] =JdT- j , (A.5) 

s 

and the solution of the Euler-Lagrange equation between (x, s) and {y, t) is 

y — ^ 

2:classical(r) = T^" - s) + X. (A. 6) 
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Inserting it into the action, we recover directly the free propagator (IA.3I 1. by imposing 
the normalization 

dyG{y,t\x,s) = l. (A.7) 



Appendix A. 2. Free particle with a wall 

We impose that some region of space cannot be visited by the particle, for example, we 
want to keep our particle in the positive region a; > 0. To do this, we put a potential 

{+00 if a; < , 
II X > . 

The solution of the associated Schrodinger equation is 

[2 

ipk{x) = \ -sm{kx), (A. 9) 

V TT 

with < A; < +oo, for x > 0, and ipk{x) = for a; < 0. Now the normalization constant 
is found directly from what we impose on the free particle solution (the only freedom is 
the phase of the solution, who is unphysical). The associated energy is 
If the particle is constrained to stay below some value M, then 

Mx) = y^sin(fc(M - x)), (A.IO) 

with < k < +00, for X < M and ipk{x) = for x > M. 

Let us compute the propagator when the wall is in x = 0: 

+ 00 

Go{y,t\x,s)= j dk ^sm{ky)sm{kx)e-^'-'^''^/^ (A.ll) 



1 



e 2(t-s) -e 2('--=) . (A. 12) 



v/27r(t-s) 

The asymptotic form, when t ^ s, of this propagator is just 

lim 6*0(2/; ^l^^; s) = 6{y — x) — 6{y + x) . (A. 13) 



t- 

Appendix A. 3. Particle in a box 
The Hamiltonian is 

1 d^ 

Ho,M = -2d^ + Vo,m{x), (A. 14) 

with the box potential 

-00 if X < , 

Vo,Mix) = { ifO<x<M, (A.15) 

+00 if X > M . 
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Go,M(y,ik,.) = 2^ sm xj exp ( -(^ - 1 .(A.17) 



The solution is indexed by a positive integer n— 1,2,... 

/ 2 /nTT \ 

with energy £■„ = r?'K^ j (2M^) and the propagator is 

n=l 

Appendix A. 4- Linear potential and a wall: Airy potential 

The name 'Airy potential' is used for conveniency, and like in the main text, we use 
in this subsection the index 'Airy' for all quantity referred to the quantum mechanical 
problem of a particle in the potential 

, ^ I +00 if —00 < X <0 , ^ ^ 

w^)=r , (A.18) 

\ AX llO<X< +00 , 

where A is a real positive parameter. We want to solve the stationary Schrodinger 
equation 

HAhy1pn{x) = ~2~dx^^^'' ^Airy(2;) 1pn{x) = E^tpnix), (A.19) 

Because of the confining potential, we already know that the energy spectrum will be 
discrete, thus we use an integer n index. For x > 0, one can re-write this equation 

-2x(^x- Mx) = 0. (A.20) 

Putting il>n{x) — x{x — En/X), then xiu) — </^((2A)^/^y), one obtains 

^"{u) -u!f{u) = 0, (A.21) 

which is solved by the Airy functions (Ai(x) and Bi(,T)). That explains the name 'Airy' 
used since the beginning. We eliminate Bi(a;) because the wave function need to tend to 
zero when x — > 00. The solutions are of the form '0„(a;) = (l/iir)Ai((2A)^/^(a; — En/X) 
with K a normalization constant. The energy levels are determined by the condition 
that the wave function vanishes in x = 0. They are then relied to the zeros — of the 
Airy function Ai(x) on the negative axis («„ > 0): 

Aif-^Ej=0 ^ = a, 2^1/3^2/3 (^=1^2,...). (A.22) 

The normalized wave functions are: 

Ai[(2A)^/3:r-n„] 

^/J-Ai^[(2A)V3^-«„]dy 

Ai pX)y^x-an] ^ ^^ 23) 

v/(2A)-V3(Ai'[-«„])^' 
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where from the first to the second line, we have just modified the expression of the 
normalization constant, using the property that Ai^(z)dz = (Ai'(— «„)) . Let us 
compute the propagator of the particle in such a potential. 

GAiry{y,t\x,s) = (i/|e-(*-^)^^-|x) (A.24) 



^'^n{y)^n{x) e" 



n=l 
oo 



^ ^ Ai [i2Xy/'y - a,,] Ai [(2A)V^x - a„] ^_(,_,),„2-v3,./3 
n=i y/(2A)-V3 (Ai'[-a„])' ^(2A)-V3 (Ai'[-aJ)' 
The general form of the Airy propagator, for the potential VAiry(a;) = Vo(a;) + Ax, is 
GM^,s) = p,)V3f;A.[(2A)V3,-.„] A.[(2A)V^.-.„] ^_„„,._.,-.,3,,. 

n=i (Ai'[-a„]) 

(A.25) 

Appendix A. 5. Linear potential and asymptotic form of the Airy propagator 

For the purpose of the main text, we have to obtain the asymptotic form of the Airy 
propagator in Eq. (1A.25I ) when |t — s| — > 0. Here, we consider the simple case of the 
linear potential (constant force), we compute exactly the propagator and then we take 
the limit, showing that this is equivalent to taking the limit of the Airy problem. First, 
the action reads, between (x, s) and {y,t) along the path x{x): 

S[x{r)] = fdrl^- (^) ' + Ax(r) I . (A.26) 



Because the Lagrangian is a quadratic form in its variables (x(r), x(r)), we know that 
the propagator is of the form 

Gun{y, t\x,s)= F{t - s) e-^[^'-i— , (A.27) 

where XdassicaK''") is the classical path between (x, s) and {y, t). Without loss of generality 
because of time translation invariance, let us take s = 0. The classical path is 

a;classical(r) = ^At^ + S^(y - x) - ^At^j y + X. (A.28) 

The action, computed along the classical path, is then 

^[a;ciassicai(r)] = -^XH' + ^^^^ + ^At(x + y) (A.29) 

Because we are dealing directly with probabilities, we can compute the function F{t) 
by the normalization condition 

r" + 00 

dy GiUy,t\x,0) = l. (A.30) 

— oo 

Thus we obtain 

^ / , ^ 1 f Kt-sf (y-xf Xit-s), A , 
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For an infinitesimal diff'erence in time, i.e. for t = s + e, we have 

„ , I , 1 / A^e^ iy-x)"^ Ae, ,\ ^. 

Gvm{y, s + e|x, s) = — = exp —{x + y}]. (A.32) 

V^Tre V o Ze Z j 

Here, let's remark that A can have a dependence in the time diff'erence e (which can 
have any value at this stage). That's not introducing any dangerous time-dependence: 
we give a linear potential A(e)x for the particle at point x, with a slope given by any 
function of e, and we ask what is the probability of finding at certain point y at time 
t = s + e, given that the particle started at point x at time s. However, we see in the 
previous formula that the expression of the propagator becomes ill-defined when e — >^ 
as soon as A(e) oc e^*^^+"\ with a > 0. 

Now let us take A = p/e, with p a fixed positive number. We have at leading order 
in e < 1 

GUv. s + e|x, .) = -i= exp ("^^^ - + 1/)) + 0{t^l^). (A.33) 

The limit e — > can be taken in formula ( 1A.32I1 . and one obtains 

limG'iin(?/, s + e|x, s) = 5{y - x) e'^^. (A.34) 

e— >0 

The difference between the linear propagator and the Airy propagator is that there 
is in the latter a wall in x = 0. But if we have to take the Airy propagator between two 
infinitesimally close points, the particle, placed in x > does not have the time to feel 
the wall. Indeed, for a time e, the Brownian particle will explore a space region of order 
e^/^. So, as close as we are from the wall, we can choose an e small enough not to feel 
the wall. In this limit, the Airy propagator is equivalent to the linear propagator, that 
is, for t = s + e, 

limG'Airy(2/, s + e|x, s) = \imGrm{y, s + e\x, s), (A. 35) 

for the same expression for the linear part of the potential. With the potential 

VAiry(a;) = Vo(x) + ^x, we thus have 

limGAiry(2/, s + e|x, s) = 6{y - x) e"^'^. (A. 36) 

Explicitely, this yields the identity 

Ai[(f)^/^-a„] Ai[(^)^/^-«. 



lim 



/9r,\ ^/•^ °° Al — V — a„ Al — X — an 

I ^ ) V L^^^ J X-^.2-^r.,^r.,v^ ^ _ ^ 

V ^ / n=i (Ai'[-a„]) 



-px 



(A.37) 



Appendix B. Another calculation of the normalization constant 



In this paragraph, we compute the normalization constant Z^^ involved in Section [3] 
using the Random Acceleration Process, which is described by the stochastic equation 
of motion x{t) = ({t) where ( is Gaussian white noise. From the joint distribution of 
heights P[{/i'^}], one can determine the marginal distribution of one height, say h'^{xo), 
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for a particular Xq. The joint distribution of the heights is Gaussian and therefore 
the marginal distribution of a single height must be Gaussian (by virtue of the central 
limit theorem). One simple way to compute is to consider the marginal distribution 
of h^{nL) : let us call p{v) = Prob(/i'^(KL) = v) this probability. p{v) is a centered 
Gaussian with a variance cx^ which we now determine. 

Integrating the joint distribution P[{/i'^}] over all paths (bridges, i.e. with pbc) 
such that h'^{0) = h'^{L) = u for u e] — oo, oo[, and keeping h^{nL) = v fixed, one 
recovers p{y) : 

/l«(L) = M 

du / m''(a;)P[{/i"}]5(/i''(«L) -^;) = ^=e-"'/(2'^'\ (B.l) 

-oo J V27r(T^ 

?t'=(0)=ti 

The path integral in this last expression (IB. Ill reads explicitly: 

h''{L)=u 



Z.. J- 



oo 



fe«(0)=M 

As h'^lnL) = t> is fixed, and because the process is Markovian, one can cut the whole 
path integral into two independent blocks : one for the time interval [0, kL] and one for 
the interval [kL,L]. This yields 

1 P + CO P ^ 

p{y) = ^ dw / m"(x)e-2/.Ld-(9.^^)^ 

Z^ J —oo J 

h'^ {kL)=v 

X j Vh%x) e-y^'' 5 (^J^ h%x)dx^ (B.2) 

h'^(0)=u 

The block for [nL, L] is simply the propagator of the free Brownian motion, but the path 
integral for [0, kL] contains an additional delta function, which constraint the Brownian 
trajectories to have a null area. In the latter, let us write r(x) = h'^{x')dx' , so that 
h'^{x) = (dr(x)/dx) = r(x), and 

h'^{KL)=v 



j Vh^{x) e-^^o''^-(a^h^)' S h^(x)dx 



/l'=(0)=M 

r{KL)=v 

Vr{x) e-Uo''^dx(&r(x))2 ^ ^^(^^^^ _ ^g_3^ 

r(0)=u 

We recognize the propagator of the Random Acceleration Process (RAP) between 
position r(0) = and speed f(0) = u and position t^kL) = and speed r^nL) = v. 
The propagator Grap is the probability of finding the randomly accelerated particle at 
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point X2 with speed V2 at time t = T, knowing that it was in point Xi with speed Vi at 
time t = 0, and is given by the formula 



G'RAp(a;2, V2, T\xi, vi, 0) = exp (^-^{x2 - xi - viT){x2 - xi - V2T) - ^{v2 - vi) 

(B.4) 

Hence 

+00 

p{v) = ^ / du{Gfree{u,L\v,i^L) X Grap{0,v,kL\0,u,0)} . (B.5) 

— 00 

Using formula f lB.4ll . one is reduced to compute a gaussian integral, and one finds 

/X 3 6 ,,.2 , , 

We can finally identify the second moment of the marginal distribution 

= {{h-{KL)f) = (B.7) 

and the normalization constant 

1 



KL-\/27rL' 

which is what we found in the main text using a different method (cf. Eq. (l26l) ). 



(Bi 
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